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1. Introduction 

The density functional theory (DFT) [H |2] is one of successful frameworks in the theory 
of electron systems. For a present standard scheme of this theory, the universal energy 
density functional F[n] defined by the constrained minimization method [H HI [5l [6] is 
used as a key quantity. Lieb has analyzed F[n], which is often cited as the Levy-Lieb 
functional. [7] This functional connects a Fermion density n(r) to a variational energy 
of an electron wave function ^> showing the same density. Correspondence from n(r) to 
this minimizing state vector |\E') is shown to exist, [7] which is a basic principle in the 
density functional theory. 

A continuous modification of the Coulomb operator appearing in the definition 
of F[n] is one of relevant techniques to analyze this functional. [H [9l O |6] The idea 
is similar to the ordinal perturbation theory of electron systems. |10l [TT] To have an 
exact expression of the so-called exchange-correlation functional, people fixed the charge 
density and considered continuous reduction of the coupling constant by multiplying a 
factor A in the energy density functional. [H [9] 

We re-analyze this modified functional, which is called F\[n] in this article, to 
consider phase transition points appearing in a generalized Kohn-Sham scheme. [I2] We 
have two main purposes for this analysis. The first purpose is to study a parameter 
differentiability of this A- modified F[n]. Fixing a density at an iV-representable density 
in Iiv,[7] the parameter derivative of F\[n] is shown to be well-defined. Namely, existence 
of state vectors giving Dini's derivatives is shown. As the second purpose, a condition 
on phase transitions around the f-representable density is derived from the theorems 
on -F\[n]. Existence of level-crossing points, i.e. phase transition points, is known. We 
address existence of an e-vicinity around the true solution of the many-electron system. 
In this finite region, no level crossing point is found provided that the density of the 
Coulomb system is also f-representable with a model given by -FaM- The definition of 
the model is given in section |6l Relevance of this £- vicinity for a general Kohn-Sham 
scheme is discussed in section El 

2. The setup of the problem 

We consider a static state of a material. The number of electrons is fixed to be a 
finite integer. We apply the Born-Oppenheimer approximation (BOA), in which motion 
of nuclei is separated from the motion of electrons. To find a stable state, we consider 
a classical state of nuclei and fix the coordinates of nuclei, whose number is M. Motion 
of the electron system is determined for this classical configuration of nuclei, which give 
a static electric potential for electrons. 

Let |0) be the electron vacuum. The electronic state in an external scalar potential 
is symbolically written by a state vector \^). This electron state is described by Fermion 
field operators, and ipa{^)) satisfying the canonical anti-commutation relations. 

K(r),^.'(rO], =5(r-r')5.,... 
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Considering a position basis | {rj, cTj}) = YliLi V'ii(ri)|0), we have a wave function as an 
inner product between this position vector and |^). 

*({r,,(j,}) = ({r,,(j,}|*). (1) 

The number operator n(r) probing existence of an electron at a point r is defined 



as 

n{ 



r)^J2'^Ur)Mr). (2) 



The kinetic energy operator T for electrons is assumed to be, 

with the electron mass m. 

Owing to BOA, we have an external scalar potential Vext(r) given by the charge of 
nuclei. Using position vectors R/ (/ = 1, • • • , M) of fixed nuclei, the potential Vext(r) is 
given as. 



7=1 



^^y^I^ (4) 

T 1 I -t 



Here, the charge of the /-th nucleus is Zie. The potential term for the Hamiltonian of 
electrons is given by the next operator. 



= j dVvext(r)n(r). (5) 

In the following discussion, we omit a constant Coulomb energy coming from the ion-ion 
interaction. 



Rf — R . , 

(i,j) I ' 

although this term is important for the charge neutrality condition. Hereafter, the 
symbol (/, J) denotes a pair of two different integers ranging from 1 to a finite integer 
(here it is M), and the summation with respect to these pairs are written as j)- 

If we have more than two electrons, the electron-electron interaction always takes 
place. In a static state, the inter-electron interaction is described by the next operator. 



a, a' 



^^Jd'r d'r' y—^^ : n(r)n(r') : . (7) 

The symbol : O : for an operator O denotes the normal ordering of the field operators. 
In : O :, the order of ijj and ijj^ is reorganized so that the creation operators come to the 
left of the annihilation operators by interchanging operators. If interchange of two field 
operators is made m times, the sign of (—1)™ is multiplied to the reordered operator. 
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We consider an isolated electron system having > 2 in an equilibrium. 
Competition between the Coulomb interaction V^e and other single-particle parts, 
T + V^xt, causes various phase transitions in the electron systems. We search for the 
minimum energy allowed for the electron system. Thus any state vector, which has an 
anti-symmetric property for the spin-1/2 Fermion system, has to be normalizable, and 
has to have a finite kinetic energy. The class of wavefunctions for the allowed state 
vectors is H^{B?^). 

For any state vector of an N electron state with > 2, we have = 
iljai{j')'ipaij)\^) 1 and the quantum state has a positive semidefinite norm, : 
n(r)n(r') : |^) = (\1>'|\I^') > 0. If we assume that (^iKel^) = 0, positivity of the 
Coulomb kernel requires that Vr and Vr', (\E'| : h{Y)h{r') : = 0. But, then we have 
zero of the integral 

j d^rdh\m\ : n(r)n(r') : |^) = ^) = q, 

which yields = or = 1. 

The expectation value of Ke by \E' G H^{B?^) is known to be finite. [7j Actually, 

I- ^ 2 

(^iKel*) = E / \{dri-^—-^*{{rk,crkm{{ry,ak'}) 



^5: J 


r ^ 

1=1 








'JW' 


„ A 


T 


+ \^drie\\ 


= Tij and 9{x) - 



N 



r.,|)]|^(K,(7,})|^J><oo. (8) 



A^ > 2, we have, 

< (^iKcl^) < oo. (9) 



3. The charge density being the primary order parameter 

The electron system is characterized by the electron charge density n^,(r) = (\E'|n(r)|\l/). 
This physical quantity is observable by X-ray diffraction measurements. Since the wave 
function \l/({rj, cTj}) is given as a function in the space of integrable functions, and 
since we assume that the kinetic energy is finite for r;,,j,(r)^/^ and Vn,i,(r)^/^ is in L?'. 
Then, n^,(r) is known to be also in -^^'^.[7] The space L^/^ + is a dual of the + 
space. Since the Coulomb potential is in L^/^ + fext(i")n^(r) is integrable. 
The Hamiltonian of this electron system is given as 

^ = f + v;e + Kxf (10) 
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The total energy E\^,Vcy±\ of this electron system is given by, 

E[^,t;ext] = (^1^1^) = (^|T + V;e|^) + j rfW(r)nv,(r). (11) 
The lowest steady state is given by minimizing this energy in the space of wavefunctions. 



EQ[vey±] = minE[^,fext] = min 



(^|f + V;e|^)+ / rfVt;ext(r)n*(r 



.(12) 



A definition of the order parameter is given by a derivative of the total energy with 
respect to the external field. For the electron system, this external degrees of freedom 
is given by fext(i')- Let's consider a stable electronic state |\E') without degeneracy. This 
state has an energy £^obext] = (^l-f^l^)- If an infinitesimal variation of fext(r) is given 
in L^/^ + the wavefunction minimizing E[i!,Vext + ^v^xt] exists. This minimizing 
wavefunction may be given as 

1^ + 5^) = e^^(l^) + |5^)). 

A complex phase factor e*^ represents a gauge degree of freedom. We also use a notation. 



6H = J d-'r6vext{r)n{r). (13) 

When ||5fext|| = 5 — )■ 0, the norm of |5\E') goes to zero, and the state |5\E') is orthogonal 
to Then, we have a derivative of i^o^cxt]- 

Eo[Vext + SVcxt] - Eo[Vcxt] 



d'r5vUv)n^ir) = / d^r^-^^^dv^r)- (14) 



This derivation follows a proof of the force theorem by Parr. [13] 

In a steady state, since the state has to be stable against any perturbation, 
the functional derivative should exist. The coefficients of the functional derivative 

- — = n^(r) is the primary order parameter of the electron system. If we have a 
dVext(r) 

first-order phase transition by introducing 5^6x1(1"), a level crossing in the lowest energy 
state |\E') occurs. Then, a jump in directed derivatives could be found. In such a special 
point, the derivative becomes ill-defined. However, this jump is detected as a jump in 
the charge density. In this sense, the electron charge density should be the primary 
order parameter. 

In sectional we will define another phase transition without a jump in n{r). This 
transition can happen in electron systems owing to internal degrees of freedom, e.g. the 
electron spin. The transition without change in n(r) can happen much frequently, when 
the interaction strength or the form of the inter-particle interaction are modified. Once 
the interaction strength is shifted from that of the Coulomb interaction, the resulting 
Hamiltonian describes a model system. The model can be identical to the Kohn-Sham 
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model. [2] Thus, for the analysis of the Kohn-Sham scheme, it is important to analyze 
this second-type phase transitions in an abstract model space. 



4. The density functional theory as a Landau theory 

To define the universal energy density functional, a Sobolev space Xn of functions in 
M'^ was introduced in the density functional theory. For a function n{Y) in the set X^v, 
nir) > and n^/^(r) is square integrable. Its gradient Vn^/^(r) is also square integrable. 
In addition, nfr) satisfies, 

d^rnir) = N. (15) 

We use the theorem 3.3 of Ref. [Tj stating that, if n{r) G Tn, F[n] given by the next 
definition exists. 

F\n] = min(^|r + Kel^)- (16) 

The symbol ^ — )■ n represents that a minimizing state |\E') is searched with a constraint 
(\E'|n(r)|\I') = n{r). Then, the existence of the minimizing is relevant. 
We have the next constrained minimization process. 



-^obext] = minminE[^,t;ext] = min 



F[n] + j 



d\vc^t{r)n{r) 



(17) 



Here, F[n\ behaves as a free energy of the electron system. 

If we have a well-defined description of the free energy of the system as a functional 
of the primary order parameter, we have a complete expression of the Landau free 
energy. The density functional theory actually gives an example. In an A^-representable 
form of the energy density functional, the Landau free energy of the system is given as 
a summation of the so-called universal energy density functional and the energy of the 

external scalar poteutial. The latter. / d'r.„(r)„.(r). is Imear m the order parameter. 
5. The A modified functional 

Let A be a real parameter in [0, 1]. We now consider a reduced energy density functional 
F\[n] defined by, 

Fx[n] = min (^'|f + XZ,\^'). (18) 

Existence of the minimizing \E'' is given also by the theorem 3.3 of Ref. [7J. We would 
like to address a next statement. 

Theorem 1 //O < A < 1, and ifn{r) G Xjy, Fx[n\ is a monotone increasing continuous 
function of X. F\[n] is concave as a function of X. 
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First, choose < Ai < A2 < 1. Assume that -^aJ^] > F\^[n]. Choose a minimizing 
state |\I/o) — 7- n of the expectation value (\1/'|T + A2Kc|^')- Then, we have, 

FaJti] = min(*'|f + AiV;e|^') 

> min(^'|f + A2V;e|^') = (^o|T + A2Kc|^o) 

= (*o|T + AiKel^o) + (A2 - Ai) X (^olKel^o) 

> (^olT + AiKcl^o). 

This inequahty contradicts the definition of Fxjn]. Thus, Fx[n\ is a monotone 
increasing function of A. Next, assume that F\[n] is not continuous when A = 
Aq > 0. This is equivalent to a statement that > 0, V5 > 0, 3A > 0, 
s.t. {|A — Ao| < 5 & |FA[n] — -FaoMI > For simplicity, let us further assume that 
< A — Ao < 5 and then Fx[n] — FaqM > ^- If "^ij) ^ '^Ni for any minimizing state 
vector |\E'o) of -FaqM; the wave function of \E'o — )■ n(r) is in if^, and 

< (^olKcl^o) < 3Co < 00. (19) 

If we let 5 = e/ (2Co), we have, 

min(^'|f + AKcl^') 

> min(^'|f + AoKel^') 

= (^o|T + AoKe|^o)+£ 

> (^o|T + AKel^o) + (Ao - A)Co + e 

> (^o|r + AKe|^o)-5Co+£ 

= ($o|T + AKe|^o) + |. (20) 

This inequality contradicts the definition of the minimum. 

Consider < Ai < A2 < 1, < ^ < 1, and A^ = ^Ai + (1 - 0^2- We call a state 
^A^ which minimizes (\E''|T + A^t^eel^')- Then, we have 

i^5Ai + (l-OA2 M 
= e(^A,|T + AiKel^A,) + (1 - 0(^A,|T + A2Ke|*A,) 

>eFA,M + (l-O^A.M. (21) 

This inequality ensures concavity of -Fa[^] as a function of A. ■ 

Following knowledge on the monotone increasing continuous functions and the 
convex (concave) functions, we immediately obtain results on derivatives and an integral 
of the derivative. |14] Let's define Dini's derivatives of -FaI'^], 

D-iX) = hm sup ^^^^N -i^AM ^ ^22) 

y-^^- y<h<0 h 

D^{X) = hm inf ^A+^N -i^AW^ (23) 

y^O+0<h<y h 

When A = 1, we formally define -D+(A) by considering Fx[n] for A > 1. In practical 
simulation, however, -D+(A) is not required at A = 1. 
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Corollary 2 For n{j) G X^r, the monotone increasing concave function F\[n] of X has 
a directed derivative at VA G (0, 1] and -D^(A) > -D+(A) 

When D~{X) = D+(A), we have the derivative, —Fx[n], of the monotone increasing 

uA 

function of -FaM- The differentiability of a monotone increasing function is given by H. 
Lebesgue. Besides, the next statement holds. 

Corollary 3 //O < A < 1, for n{r) G X/v, Fx[n] is differentiahle a.e. The derivative 
— FA[n] is integrable in [0,1]. 

Next, consider the case with > 2. A set of state vectors |\E') reproducing n(r) 
is denoted as Q„. Fix A G (0,1]. We consider a subset of |\E') G Qn minimizing 
(^|f + AKcl^) and call it P„(A). For any |^) G P„(A), < (^l^eel^) < oo. Thus 
we have a finite range including values of (\E'|\4e|^) for |^) G P„(A). The maximum 
of this range is given by a state vector in Vn{X), and also the minimum is given by 
another vector. They might be different with each other. Thus they are denoted as 
l^^x) ^ 'Pn{X) and |\E'^) G 'Pn(A). The very definitions of 'Pn(A) and these vectors ensure 
that 

Fx[n] = (vl/-|f + AKel^A ) = (^A 1^ + >^tcM), (24) 

and that 

i^tlVeel^t) < (^A iKel^A)- (25) 

Now we prove existence of a vector giving the directed derivative of -Fa[^]- 
Lemma 4 For A G (0, 1], we have D+{X) = (^+|V;e|^^) and D~{X) = (^;^|V;e|^A )■ 
For < V(5 < A, we have 0<X' = X- 6<X and 

Fy[n] < (^;^|f + AVec|^A) = (^aI^ + A^ccI^a) -'^(^Al^'eel^A) 

< (*+|f + XtelH) - <^(^A l^eel^l) = (^^ |T + AVeel^;^). (26) 



Thus, 



n-/^^ r Fx[n] - Fx^h[n 

D (A) = lim sup 



y~*o+o<h<y h 

> lim sup — 

y^o+ o</i<j/ h 

= {^xlV^el^x). (27) 

Similarly, we have, 

I^+(A) < i^tlVeel^D- (28) 

These inequalities together with Eq. fl2^ give just another proof of Corollary [21 
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Let's assume that D~{X) > (\I/^|Vee|\I/^). Then, we have 

F,[n] - F,^,[n] - h{<f-\te\^-) 
< hm sup — 

0<h<y h 

= lim sup (^Ilr+(A-ft)V;.|^,-)-FWn], ^9) 

y^0+ o^h<y h 

Independently, we have a next inequahty by the definition of Fx-h[n]- 

(^-|f + (A - h)Vee\'^l) > Fx-h[n], (30) 
which yields, Wh G [0,y], 

f{h) = {^-,\f + (A - /i)v;el^A ) - ^A-^N > 0. (31) 

By Eq. ([21, /(O) = 0. 

We show that f{h) /his also a monotone increasing function. Actually, if we assume 

that < 3hQ < y, and that > li^^ we have a next inequality. 

ho y 

-h^F,[n] + ^F,_y[n]>F,_,M- 
This contradicts to the concavity of Fx[n] given by Eq. f l2T]) . Thus, Eq. fl29|) tells that 

< lim in\f+(X- v)K.K)-F.-.ln]^ ,32, 

y^o+ y 

We call a state vector, which is in 'P„(A — y) and maximizes the expectation value of 
Vee, l^^x-y)- Then the above expression yields, 

y—>-0+ y L " ^ 
> inlVeel'ifx). (33) 

We note that limj,^o+ I^a-j/) exists. The continuity of Fx[n] ensures that 
limj^^o+ l^^x-y) ^ 'Pn{X). The inequality, Eq. f l5^ . contradicts to the definition of l^'^). 
Thus we conclude that D-{X) = (^"iKel^A)- Similarly, D+(A) = (^+|V;e|^A)- ■ 

For a fixed n{r), we now analyze the number of discontinuous points, where 
D{X) = D-{X) - D+{X) > 0. Lemma H tells that, when D{X) > 0, P„(A) has muhiple 
elements which are distinguished by difference in the expectation value of Vee, and 
also in the expectation value of T. At this discontinuous point, we have a change in 
the minimizing state from I^E'^) to |^;^) by reducing A. Therefore we have the next 
definition. 

Definition 5 A point where -D(A) > is called a transition point in the model space. 

We should note that a priori I^P^) and I^E'^) are defined independently with each 
other. Lemma m implies that, for a point with -D(A) = 0, we can identify I^E'^) with I^E'^) 
and replace one with the other. Thus the point with D{X) = is not a level crossing 
point. 
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Theorem 6 i) In [0,1], we have at most a finite number of transition points, where 
F\\n\ has a finite discontinuity, which is greater than an arbitrary small number e. 
ii) At an accumulation point of the discontinuous points D{\i) (i = 1, 2, ■ ■ ■ , ooj with 
Xi G (0, 1], and Aj < Aj+i, we have -D(Aj) — )■ for i — )■ oo. 

i) For any finite number e > 0, we have a set A{e) of discontinuous points, at which 
D{X) > e. Assume that the number of elements of A{e) is more than the countable 
infinite. We can select a countable infinite subset of A{e) and name it A.{e). Points in 
A{e) are in [0, 1] and are to be ordered as Aj < Aj+i. This is due to the selection axiom. 
Let's number the points in A{e) as Aj (i = 1, 2, ■ ■ ■). We have. 



Since the set A{e) is infinite, we have an integer J > D^{\i)/e for which the derivative 
Dj^{\j) becomes strictly negative. This result contradicts to the increasing property of 
^^[n]. Thus the set A{e) has to be a finite set. 

ii) We consider infj D{Xi) for i = 1,2, ■ ■ ■. Owing to i), a finite non-zero infimum is 
denied. Now consider D = lim7v-s.oo supj>^ D(Ai). Assume that D > 0. Then, we can 
find an integer li, for which -D(AiJ > D. But, we can also find another integer I2 > h, 
for which D{Xi^) > D, because of the definition of D. Then we have an infinite series 
of A;, (i = 1, 2, ■ ■ ■) with D{Xi-) > D > 0, which contradicts to i). Thus, D = 0. M 

Lemma 7 For n{r) G In, th^ function ^-^^M of X is Lipschitz continuous in [0, 1] 
and we have. 



For V\[' G H^, we have a finite R satisfying (\E'|iyee|^) < R- Choose R such 
that -D+(0) < R. Since Fx[n] is continuous monotone increasing and concave, for 
< VAi < VA2 < 1, we have. 



< Fa.M-Fa.N < D+(Ai)(A2-Ai) < D+(0)(A2-Ai) < i?(A2-Ai).(35) 
Thus, we have Eq. (1341) . ■ 



6. The Kohn-Sham minimization scheme 

To discuss relevance of our theorems for the discussion of ground-state phase transitions, 
we introduce a Kohn-Sham minimization scheme. [T6| [T2] For an external potential 
fext(r), we have a ground state of the Coulomb system, I^E^gs); which gives a ground-state 
density, ?T,Gs(r). First, we have a next equality. 



D+(A,) < D-(A,) < Z^+(A,-i) < /^-(A,„i), 

D+(Aj) < D-{Xj)-E< D-{Xj^i)-2e < D-{Xi) -je. 




(34) 



= (^Gs|T + V;e|^Gs)+ / 
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mm 



mm 



mm 



mm 



F[n] + j (i^rfext(r)n(r) 



mill (^'|T|*') + Fx=i[n] - Fa=oN + / c^^rfext(r)n(r) 



mill <j (^'|r|^') + / d\—F^[n^,] + / dhv,^,{v)n^,{v 



d 



d\^Fx[n^,>] + / d^rv^^t{j)nm,{r) 



minGc.^^J^^']. 



Here, we used a notation for the charge density given by \^) as, 

n,i,(r) = (^|n(r)|^). 
We can show two conditions on the minimizing \1/ of GctJcxt [^] • Note that, 

Go,...J^] = (^|r|^) + Fa=iK] - Fa=oK] + j dhv,^,{r)n^{v) 
Then we have. 



(36) 

(37) 
(38) 



(^'|T|^') + Fx=i[n^] - FA=o[nv,] + / d'rv,^,{v)n^{Y) 
= Fa=i + y dVwext(r)^*(r) 

>i^A=i[riGs] + j dhve.tir)nGs{r). (39) 

We see that equahties in Eq. ([32]) are satisfied, 

(i) if \l/ is identical to a states which minimizes (\I/'|T|\[^'), with the constraint that 

(ii) and if is identical to the true ground-state charge density r;,Gs(r). 

Thus, the minimizing process of G'o,t,^„J\l'] gives us a state that reproduces nGs(r) and 
minimizes (\1/|T|\E^). Two fundamental statements are addressed here. 

(i) A minimizing state \E' of Go,i,ext [^] is the state which is searched in the constrained 
minimization of -Fa=o['^gs]- Thus is the state motivated to be searched in the 
Kohn-Sham scheme. 

(ii) We do not need to have a secular equation to define the minimization process of 
Go.iiext [^] for our discussion. 

The definition of Go.j^ext [^] suggests us that we have a plenty of models for electron 
systems. Actually, we have another equality. 

Eo[ve^t] = mmGA,^,ext[^], (40) 
G'w[^] = (^l^ + AKc|^)+ f' dX'^Fy[n^]+ /"dW(r)n*(r).(41) 
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The wave- function funcional, G^.^extl^]) determines I^I/a), which may appear as a 
minimizing state in the definition of -Fa[^gs]- Therefore, when we fix fext(r), 
minimization of G'A,t;cxt[^] with A G [0,1] produces a set of states {|\I/a)} and thus 
the function -FA[nGs] of A as 

i^A^Gs] = (^aIT + AKcI^a). (42) 

We may introduce the Hartree term to formulate a much famihar form in the density 
functional theory. 

Go,.... m 

+ [ i dx4r min (^\f + XVJ"^) --[ d^rd^r'—^ -n^(r)n^(r')\ 

Jq { c/A *'^n* 2 J |r-r'| J 

= (^ITI^) + - f d^rd^r'—^ -nq,(r)nq,(r') + f (iVt;ext(r)n^(r) 

2y |r-r'| J 

1 r 

+ / dX- / d^rd^r'- :{{^x\ ■ n(r)n(r') : |^a) - n^(r)?2^(r')} 

Jo 2 J |r-r'| 

= (^ITI^) + - / d^rd^r'—^ -nMn^ir') + / d^rv^^ti^c^n^ir) 

2 J |r-r'| J 

+ E^,[n^]. (43) 
Similarly, we have 

= (^|f + AKecl^) + j d^rd\'^^^—^^n^{r)n^{r') + J d^rv^^t{r)n^{r) 

+ / dX'- d^rd^r'- {("^ ■ H^)n(r') l"^ y) - n^(r)n^(r')} 

Jx 2 J |r-r'| 

= (*|f + AKel^) + y ci\rfV^-^--^n^(r)r2.i,(r') + y rfVi;ext(r)n^(r) 

+ ^xc,A[r^vi.]. (44) 

We call Gx,vcy,t ^^e A-parametrized model energy functional. All of these models can 
determine i^of^ext] as shown by Eq. f HOj) via a determination of |^a)- The charge density 
of |\E'a) satisfies, n^^(r) = 'n,Gs(r). When < A < 1, since a reduced- interaction AV^e 
appears in the definition of Gx,va^^[^], its minimizing state, l^^x), is represented by a 
summation of Slater determinants. 



7. Existence of e vicinity 

From now on, we consider possible level crossings in -FAfncs] given in Eq. ( l42i) . A level 
crossing occurs, when two or more states appear as minimizing states |\I'a) of T + AV^g 
and when a jump in Dini's derivatives, -D(A) > 0, happens. 
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Existence of the crossing points is exemplified by a phase transition from the normal 
state to a ferromagnetic state recognized in the uniform electron gas system. [TOl [T9] At 
the phase boundary, we have two uniform electron gas ground states, i. e. a paramagnetic 
state without spin polarization and a partially ferromagnetic state with a finite total 
spin. In general, we need to assume that crossing points appear at some of A G [0, 1]. 
Once there appears a crossing point, we can count the number of crossing points on the 
A axis. 

First, we note that the level crossings happen at discretized points on the A axis. 
Two statements of Theorem |6] deny a possibility to have crossing points with D{X) > 
continuously or densely on the A axis. Next we need to consider a case with countable 
infinite numbers of crossing points in a finite region of A. This case causes appearance 
of accumulation points of the crossing points on the A axis. 

Thus we can choose a finite interval of A satisfying one of two possible conditions 
for -D(A): i) VA G (Aq, A;,), -D(A) = 0, or ii) one of the boundaries Xi {I = a or b) is an 
accumulation point of -D(Aj) > and limj_j.oo D(Aj) = 0. Our concern is whether the 
case ii) happens at At = 1 for a t>-representable ground state density or not. 

Consider a unique ground state |\E'gs) of an electron system in the external potential 
fextlr), which has an electron charge density ?2gs(i")- Simbolically, we write the v- 
representability of nGs(r) as nGs(i") ^ ^Ar-[Z] The ground state is supposed to be stable 
against small perturbation. Let's assume that the point of A = 1 is an accumulation 
point of the level crossings along the A axis, and that we have the case ii). 

Crossing points are labeled as Aj (i = 1, 2, ■ ■ ■ , oo) and Aj < Aj+i for the present 
discussion. At any level crossing point Aj, owing to Lemma H] there are at least two 
minimizing states, I^&a,) and I^&a.)- The expectation values (T + AjV^e) by these states 
are the same at the crossing point. However, l^^^,) and {"^x ) distinguished by (Ke) 
and (T). Thus, they are linearly independent as state vectors. We can also show that 
(\&^. I V^el^A ) decreases, when i increases, as exemplified in the proof of Theorem |6l 

In between two neighboring crossing points, we have a continuous change in 1^^^) 
and -D(A) = 0. The state, |^^.) is connected to |\I'^.^^). We may select as a 

representative state for this finite range [Aj, Aj+i]. 

Existence of a big number of level crossing points at a close vicinity of this 
accumulation point requires existence of a plenty number of nearly degenerate states, 
1^^.) with i = l,2,---,oo. Here, the state, \'^x,)y flight not be an eigen state of 
any potential problem, but minimizes just the expectation value of (T + A\4e) keeping 
the density. They are not distinguished by the primary order parameter, n(r), nor by 
difference in any external symmetry breaking observed by the external potential fext(r). 
Only internal degrees of freedom distinguish these infinite numbers of 

We should note that only to have these states does not directly mean existence 
of infinite numbers of degenerate eigen states at A = 1. This is because we have a 
possibility that all of I^E'^.) except for I^'gs) become non-eigen-states of H, but are just 
variational states. Furthermore, if continuous change in {"^x) is allowed in a finite range 
of A, only one minimizing state appears as a unique minimum for (T + AV"ee) in the 
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range. 

However, the state exists and may be used as variational states at any 

A G [0, 1]. By a simple inspection, we can see that at A = 1, the variational energies of 
(^^. |T + Vccl^A ) separated by a finite gap with each other. However, when A = 1 
is an accumulation point, we have |(^';[.|T + Veel^Aj) ~ (^X J-^ + ^cc|^Ai+i)l ~^ ^ 
i — > oo. Thus the gap ioi i^ N should be much smaller than any energy separation in 
the energy spectrum of the finite size system. 

To discuss the accumulation points further, we now analyze existence or non- 
existence of a potential v\{r) which gives a secular equation determining a state vector 
|$a)- Namely, for A = 1 — 5A with < 5A <^ 1, we formulate a quantum mechanical 
problem, whose solution satisfies ||($A|^(r)|$A) ~'"'Gs(r)lloo — 0, and |<I>a) I^gs) 
for A — 7> 1. Simbolically, we write this statement on another i>-representability of nGs(r) 
as nGs(r) G Ax,n- 

For simplicity, we consider a compact space by introducing a cube with volume 
under the periodic boundary condition. Then, we can introduce a Fourier series 
expansion for the potential v\{r), which is used as the Lagrange multiplier to fix the 
charge density of |$), as, 

^^^^) = ^^.G cxp(iG ■ r). (45) 

G 

Consider an N electron state |$). We introduce a functional Q\\^),v\^g,E : 5\\ as, 
- ($|f + (1 - 5\)V^ + / d\^;ext(r)n(r)|$) 



+ j dhvx{v) K$|n(r)|$) - nGs(r)] - E [($|$) - 1] 
= ($|f' + V;e + j dhv,^{v)n{v)\^) 

+ 1 E E E ^a,g4_g,.cg',. - 5x%, \ 1$) 

Igt^O G' fT J 

- ^A,GnGs(-G) - E - 1] . (46) 

Here, nGs(G) is the Fourier component of nGs(r), Cg,o- are electron annihilation 
operators with the spin cr, and E = E — Nvx^. By making derivatives of Q with 
respect to variables except for a parameter SX, we have next secular equations. 

[h + Hsx}\<^) = E\^), (47) 

^SX^J2J2Y1 ^A,G4._G,aCG',a - SXV,^ (48) 
Gt^O G' o- 

($1$) = 1, (49) 

E E(*I4'+g,.cg',.|*) - nGs(G). (50) 



G' a 
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At A = 1, 5A = 0, the solution of Eqs. (H7|) . ( H^ . and (ISO]) is given by the normalized 
state I^Gs) with v\^g, = for all G. For any vx{v) E L^/^ + or equivalently for any 
set of fA,G, with 6\ > 0, we have a normalized eigen state |$) of Eqs. ( H7I) and ( l49l) . So 
the construction of vx^g to meet Eq. (!50!) is the problem. If t>A,G exists, nGs(i") G ^A,Af- 

Lemma 8 When a unique ground state of H gives the charge density nQs{r), and if 
there happens accumulation of points Aj (i = 1,2, ■ ■ ■ , oo) satisfying < Aj < Aj+i < 1 
and D{Xi) > 0, nGs(r) is not in Ax^^n- We have a series of potentials fAi(r); which 
gives a set of pure states \'^Xi), and ||(<l>A-|n(r)|<l>Aj — ?^Gs(r)||oo —^0 as z —t- oo. 

To show this lemma, we classify possible conditions and deny possibilities except 
for that ?T,Gs(r) is not in ^a^.a^, which is the case 4 below. 

case 1 There exists fAi(r), and I^E^Ai) ^^"^ I^Ai) eigen states oi H + Hsx = 

f + X,Vee + j dhvxXr)n{T). 

case 2 There exist v'^_{t) and v^_{r), which are different from each other more than 

a constant. Two states, I^^a^)) are eigen states of T + XVee + J d^rvf,{r)n{r), 
respectively. 

case 3 One of I^E'a-) and I'^Xi) eigen state of T + XVe + J d^rvxi{r)h{r), and the 
other is not. 

case 4 Both of minimizing states, I^E'a ) and I^I'a not eigen states of any potential 

problem given as T + XVe + J (i'^rt>A.(r)n(r). 

We first deny the case 1. Let's assume that the solution t>A,G exists for nGs(r)- It 
means that we have a minimum of Q, which is given by an eigen state |$) of Eq. (I47p . 
If we insert I^'a) in the functional Q, we have 

Q[\^x),vx,G,E:6X] 

= (^A|T + AKe|^A)+ / c/WtWriGsW 



< mf + XVe\^) + j rfVt;e,t(r)nGs(r) 

= Qm,vx,G,E:6X]. (51) 

Thus, by the variational principle of the quantum mechanics, we see that |\E'a) is also an 
eigen state of Eq. (H7|) . If we choose A = Aj, we have two degenerate eigen states, I^^Ai) 
and I^Ai)- They are orthogonal with each other. By taking the limit i — > oo, we have 
two limiting states, I^Aoo) I^Aoo)' since these state vectors are in a Banach space 
and existence of the weak limit is ensured by jV^el^Ai) ~^ (^GslKel^Gs) owing to 
Theorem [6l Orthogonality between I^Aoo) ^^"^ I^Aoo) Iiol*is. Of course Aoo = 1- These 
states minimize Qd^E'), f a=i,G) • ^X = 0]. This means that the ground state has to be 
degenerate, which contradicts to the uniqueness of the ground state providing nGs(r)- 
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If we have a case 2, we have two independent potential v^,{v) and f^. (r), both 
of which give the same density nGs(r)- This case contradicts the Hohenberg-Kohn 
theorem [1] and thus it is denied at any i. 

If we have a case 3, we again have a difficulty. Both of and |^;^,) give the 

same variational energy for T + XVee + j d^rvx^{r)n{r). The variational principle tells 

that a state having the variational energy of the lowest eigen state gives a degenerate 
eigen state. This fact contradicts an assumption that one of these two states is not an 
eigen state of any potential problem. 

So, we conclude the case 4, which says that vx.^g does not exist, when A = 1 is 
an accumulation point. Since nGs(r) G X/v, we know existence of an iV-particle density 
matrix F, which gives the infimum of the Lieb functional. [7] F has to be that for a mixed 
state, 

r = J2c[M.){'^i\- (52) 

I 

l^^x^) is given as an eigen state of a potential problem and these states are orthogonal 
with each other. Since the limit of Aj — t- 1 is given as a pure state, when i — > oo, 
coefficients and state vectors satisfy, 

^kio, (53) 
\^t) ^ I^Gs). (54) 

Thus the second statement of the lemma holds. ■ 

The meaning of this lemma is somewhat redundant. When we have an accumulation 
point at A = 1, we have no potential series, fA,(r), whose pure ground state reproduces 
exactly nQs{r). Thus the case 4, if it is found, gives an example of an A^-representable 
density which is in An but not in A\^n for A < 1. This density is apparently not pure- 
state v-representable in any Kohn-Sham scheme. The accumulation of crossing points 
for a ground state density contradicts a picture of the ordinal Kohn-Sham scheme, 
which assumes that a unique ground state of the Coulomb problem is reproduced by a 
non-interacting system with an optimized potential. 

In this case, we should follow the present lemma to know existence of a converging 
series of quantum mechanical models. Actually, the potential series of v\.{r), whose 
ground state density is slightly different from the final solution, can be used to find a 
Cauchy sequence of nx.{r) converging to nGs(r). This converging series is found in the 
model space of the multi-reference generalization of the Kohn-Sham scheme. We may 
call the region of the model space as an e-vicinity, in which a convergence of a simulation 
is guaranteed. 

Conversely, if a density nGs(r) is pure-state f-representable in a Kohn-Sham scheme 
including a multi-reference generalization, there is no accumulation of crossing points 
at A = 1. In this case, we have a well-defined e- vicinity around the true ground state 
in the model space, in which no level crossing is found. More precisely, we have a next 
statement. If the density nGs(r) is w-representable also for A < 1, accumulation of 
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phase transition points is forbidden when A — t- 1. Finding a convergence in the density 
searched in an optimization process of model quantum systems may allow us to conclude 
no remaining level-crossing point along a line approaching in the true Coulomb system 
given by H. Therefore, we conclude the existence of an e vicinity around the Coulomb 
system in the model space, where no essential phase transition occurs in the direction 
to reduce the interaction strength by introducing A < 1, and keeping the charge density 
n(r) unchanged. 

Now we re-analyze the potential, wa.g, in a pertabative argument. Although this 
problem was treated more elegantly by Kohn, [20] we want to show a subtle problem on 
existence of fA,G- We have the limiting solution of I^E'gs) with wa,g = 0, which is an 
eigen solution of Eq. Wl\ at least as a stationary state. The construction of Q is based 
on the quantum mechanical variational principle, and the limit of Q for A — )■ 1 should 
behave regularly. In order to inspect on the w-represent ability of the normal solution 
nGs(i') for A < 1, we argue a perturbative construction method of wa,g- When the 
spectrum of H is normal, and when the excited states are written as the solution 
of Eq. (H7|) is represented by the ordinal perturbation theory as, 

1$) = ^ I I^GS) + -^^^\'^^){'^^\H5X\'^GS) + o{5\) ) . (55) 

The expectation value in Eq. (1501) gives. 



^^($|ct .CGVI$)-^GS(G) 



^ I E E ^E^:^^^'^^^' E E 4+G,.CG',.i^.) 

^ (^ilEE4'-P>-^f^"'-I^Gs)^^A,P 



G" o- 



■^-^ Eq — Ei 



(^GS I E E 4'+G,aCG',a I ^i) (^'i I ^AFee I *GS) + C.C. + 0(<5A) 

G' (T 

= E ^G.P^^P + "^^^G + c.c. + o{5\). (56) 

Here, C is a normalization constant. In general, we have a non-zero vector Bq,. Since the 
matrix Ag,p is an Hermite matrix, we have a solution fA,p of Eq. (1501) in o{5X). Thus, 
we may utilize the determination method of f a,p to analyze the Kohn-Sham method, in 
which the density nGs(r) is reproduced by another artificial model system. 

However, we find a difficulty in Eq. (1561) . The required conditions amount to the 
same number as f A,p. So, once we consider higher-order conditions for A" as independent, 
the number of the conditions is over the number of fA,p. This suggests that the solution 
could be found in a non-local potential. Even if so, the number of conditions is too large. 
If we consider all the expansions in Eq. (1561) as functions of A and f a.p, the determination 
equations form simultaneous equations of fA,p- The structure is not trivial for higher 
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order terms. Since they are non-linear determination equations, the solution is expected 
to exist only when all the expressions are derived exactly. 

Before closing this section, we summarize conditions for Lemma H] and Theorem El 
The proofs of these statements tell that the expectation value of the relevant interaction 
term has to be positive and finite, i.e., < {Vee) < oo. The convexity (or precisely the 
concavity nature of -Fa['^]) comes from linear nature of the A-modified F\[n\ with respect 
to the parameter A. The constrained minimization allows us to conclude both continuous 
nature of the functional and existence of minimizing states, which give Dini's derivative 
of the parametrized energy density functional. 

8. Summary and conclusions 

In section m we consider the A modification of F[n], while the primary order parameter 
n(r) does not change. We consider a space of models, (^A.^cxt [^]: with A G [0, 1] in section 
|6l In the model space, at a discontinuous point of Fx[n] with fixed n(r), the minimizing 
states, |\E'^) and l^^x)^ exist. Existence of discontinuous points may be detected by 
d d^ 

—Fx\n\, or ^-pt-FaM- The latter behaves as a delta function at the discontinuous 
dX dX^ 

points. We can use these functions as an indicator for discontinuous transition points 
in a set of states giving ^ — )■ n{r). 

Owing to the above discussion, we have a physical conclusion on the universal 
energy density functional. Consider a true density nGs(r), which is given by a ground 
state in an external potential Wext(r)- If the state is non degenerate (or at least finitely 
degenerated) and if it is stable, we have no accumulation point of the level crossing 
points for -Fa[^gs] at A = 1 in the searching process utilizing a potential problem with 
the reduced interaction strength. Thus, when we move away from A = 1 on the A axis, 
we have an e vicinity, where no level crossing happens in the model space, which is to be 
defined by a secular equation of the many-Fermion system with the reduced interaction. 

In a multi- reference density functional theory (MR-DFT), one of the authors 
considered general modification of the universal energy density functional. [T2l [T7] In 
addition to the Kohn-Sham description of the true density, [2] we have plenty of effective 
descriptions using model Fermion systems. The models include partially correlated 
systems, whose ground state is obtained in a multiple Slater determinants. We can 
introduce distance between two models given by the norm of the charge density and then 
the set of models with the charge distance becomes a space of models. An advantage of 
this generalization of DFT is that we can search for an optimized model in a space of 
models including non-interacting Fermion models and interacting Fermion models. 

In an optimization process in the model space, the density is searched in a 
f-representable subset. Thus, in a realization of MR-DFT, we never meet the 
accumulation of unreachable crossing points, when we see a convergence of the density 
in a searching step. The existence of an e-vicinity in the model space suggests that we 
have a physically converged model, which may continuously connect to the true electron 



Theorems on ground-state phase transitions in Kohn-Sham models 



19 



system. Even when the original Kohn-Sham model is separated from the electron system 
by crossing points, the generalized Kohn-Sham model in MR-DFT can be settled in the 
£-vicinity. 
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